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Abstract 

We study the community structure of the multi-network of commodity-specific trade relations among 
world countries over the 1992-2003 period. We compare structures across commodities and time by means 
of the normalized mutual information index (NMI). We also compare them with exogenous community 
structures induced by geographical distances and regional trade agreements. We find that commodity- 
specific community structures are very heterogeneous and much more fragmented than that characterizing 
the aggregate ITN. This shows that the aggregate properties of the ITN may result (and be very different) 
from the aggregation of very diverse commodity-specific layers of the multi network. We also show that 
commodity-specific community structures, especially those related to the chemical sector, are becoming 
more and more similar to the aggregate one. Finally, our findings suggest that geographical distance 
is much more correlated with the observed community structure than regional-trade agreements. This 
result strengthens previous findings from the empirical literature on trade. 

PACS numbers: 89.75.-k, 89.65. Gh, 87.23. Ge, 05.70.Ln, 05.40.-a 
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I. INTRODUCTION 



In the last years there was a surge of interest in the study of international-trade issues from a 
complex-network perspective [8, 9, 12, 13, 17-19, 23, 30, 36, 37]. Many contributions have explored 
the evolution over time of the topological properties of the aggregate International Trade Network 
(ITN), aka the World Trade Web (WTW), defined as the graph of total import/export relationships 
between world countries in a given year. More recently, a number of papers have instead begun to 
investigate the multi- network of trade [7, 29], where a commodity-specific approach is followed to 
unfold the aggregate ITN in many layers, each one representing import and export relationships 
between countries for a given commodity class (cf. also Refs. [21, 22] and the pioneering work of 
Paul Slater, cf. Refs. [34, 35]). 

In this paper, we explore further the topological architecture of the multi-network of interna- 
tional trade studying, for the first time, its community structure (see Ref. [15] for an overview). 
Detecting the community structure of the ITN and how it correlates with country-specific variables 
(e.g., size) and geography (e.g., distances between countries) is crucial from a international-trade 
perspective. Indeed, finding communities in the ITN means identifying clusters of countries that 
carry tightly interrelated trade linkages among them, while being relatively less interconnected 
with countries outside the cluster. To date, only two papers have been trying to explore the 
community structure of the ITN [31, 38]. However, they have only studied the aggregate ITN, 
i.e. the network obtained from total import/export relations between countries irrespective of the 
specific commodity traded. By focusing on the aggregate ITN only, one indeed neglects the fact 
that countries actually trade different lines of products and mostly employ imported goods either 
as inputs to the production process, or as consumption goods. Therefore, identifying clusters of 
countries from a multi-network perspective may be relevant to better understand what are the 
countries in the world that tend to trade the same group of products over time and, in turn, un- 
covering some stylized facts about the actual input-output and supply-demand interdependencies 
between countries. This may be relevant to predict, for example, to what extent a negative shock 
hitting a particular industry in a certain region of the world (or in a cluster) may spread and affect 
the same industry (or closely related ones) in another region of the world (or in another cluster). 

Here, we begin addressing this issue by detecting the community structure characterizing the 
commodity-specific ITN over the period 1992-2003 (T = 12 years). We employ data about 162 
countries and 97 commodities (2-digit disaggregation), to build a sequence of T multi ITNs. We 
begin by focusing on the 14 top-traded and economically relevant commodities, identifying the 
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community structure of each layer (i.e. groups of countries that mostly trade a given commodity). 
We then compare commodity-specific community structures with a number of properly-specified 
community benchmarks. These benchmarks are the community structures obtained from: (i) the 
aggregate ITN; (ii) the network of geographical closeness (i.e. the inverse of geographical distance) 
between the 162 countries; (iii) trade-related partitions obtained by detecting the community 
structure of the regional trade agreement (RTA) network. The main question we ask is whether 
(and how) commodity-specific community structures are similar to, or differ from, those detected 
in the benchmark networks. Comparisons are made using the normalized mutual information 
index (NMI), which is a measure of how close two partitions of the same set of N units are 
[10]. Understanding whether community structures detected at the commodity-specific level are 
similar to — or different from — those detected in the benchmark networks can shed further light 
on the topological architecture of the ITN. For example, comparing aggregate and commodity- 
specific community structures may tell us whether the community structure that we observe at the 
aggregate trade level can be explained by the aggregation of heterogeneous community structures 
or, conversely, trade community formation is not affected too much by the type of commodity 
traded. Similarly, comparing trade-induced communities with those obtained by the network of 
geographical closeness may help us to understand the extent to which the formation of trade 
communities is related to geographical distance (as a proxy of trade resistance factors, e.g. trade 
fees) . 

The rest of the paper is organized as follows. Section II describes the databases that we employ 
in our exercises. Section III explains the community detection method that we use in this work. 
Section IV discusses our main results. Concluding remarks are in Section V. 

II. DATA AND DEFINITIONS 

We employ bilateral trade flows data from the United Nations Commodity Trade Database (UN- 
COMTRADE; see http://comtrade.un.org/). We build a balanced panel of N — 162 countries 
for which we have commodity-specific imports and exports flows from 1992 to 2003 (T = 12 years) 
in current U.S. dollars. Trade flows are reported for C = 97 (2-digit) different commodities, 
classified according to the Harmonized System 1996 (HS1996; http://www.wcoomd.org/) [1]. 

We employ the database to build a time sequence of weighted, directed multi-networks of trade 
where the N nodes are world countries and directed links represent the value of exports of a 
given commodity in each year t = 1992, . . . ,2003 [2]. As a result, we have a time sequence of 

4 



T mult i- networks of international trade, each characterized by C layers (or links of C different 
colors). Each layer c = 1, . . . ,C represents exports between countries for commodity c and can 
be characterized by a N x N weight matrix X t c . Its generic entry x^ t corresponds to the value 
of exports of commodity c from country i to country j in year t. We consider directed networks, 
therefore in general x^ t 7^ Xj it . The aggregate weighted, directed ITN is obtained by simply 
summing up all commodity-specific layers. The entries of its weight matrices X t reads: 



For the sake of exposition, we shall focus on the most important commodity networks. Table 
I shows the ten most-traded commodities in 2003, ranked according to the total value of trade. 
Notice that they account, together, for 56% of total world trade and that the 10 most-traded 
commodities feature also the highest values of trade-value per link (i.e. ratio between total trade 
and total number of links in the commodity-specific network). In addition to these 10 trade- 
relevant commodities, we shall also focus on other 4 classes (cereals, cotton, coffee/tea and arms), 
which are less traded but more relevant in economics terms. The 14 commodities considered 
account together for 57% of world trade in 2003 [3]. 

We also employ data about regional trade agreements (RTAs) between world countries taken 
from the World Trade Organization (WTO) website [4]. We build a weighted undirected network 
with weight matrix M t = {m^} where nodes are countries and a link is weighted according to 
the number of RTAs - free, multilateral and/or bilateral - in place between the two countries 
i and j at year t [cf. also 31]. This sequence of networks may be interpreted as an indicator of 
how intense are trade agreements between countries over time, i.e. how close countries are in the 
RTA space. It is well-known from the empirical literature on trade that RTAs are an important 
determinant of trade flows [11]. 

Finally, we build a geographically-related weighted undirected network with weights s^- = d" 1 , 
where dij are the geographical distances between the most populated cities of country % and country 
j[5]. We employ the resulting matrix S = {sij} as a weighted undirected network of geographical 
closeness between countries, i.e. as the network conveying information on how close countries are 
in the geographical space. Notice that, traditionally, geographical distance between countries is 
interpreted as a proxy of all factors that impose some resistance to free trade (transport costs, 
fees, etc.). 



c 




(1) 



c=l 
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III. COMMUNITY DETECTION AND COMPARISON 



It has been observed that many real networks exhibit a concentration of links within a special 
groups of nodes called communities (or clusters or modules). Such a structural property of a 
network has also been linked to the presence of sub-modules whose nodes have some functional 
property in common. Therefore, the detection of the community structure of a given network 
could help to discover some hidden feature of its topological architecture. 

Despite the intuitive concept of community, a precise definition of what a community is rep- 
resents a challenging issue (see Ref. [15]). In this paper we adopt the well known formulation 
given in [26] : a subgraph is a community if the number of links (or, more generally, the intensity of 
interactions) among nodes in the subgraph is higher than what would be expected in an equivalent 
network with links (and intensities) placed at random. This definition implies the choice of a so- 
called "null model", i.e. a model of network to which any other network can be compared in order 
to assert the existence of any degree of modularity. The most used null model is a random network 
with the same number of nodes, the same number of links and the same degree distribution as in 
the original network, but with links among nodes randomly placed. Based on these concepts, a 
function called modularity that gives a measure of the quality of a given network partition into 
communities has been introduced in Ref. [26] . The modularity function has been further extended 
in Ref. [6] to the case of weighted directed networks as reported in the following: 



Q = —Y 



w 



Scucj (2) 



where Wij is the weight of the link between i and j, w° ut = J2j w ij an d wj 1 = J2i w ij are respectively 
the output and input strengths of nodes i and j, W = J2iJ2j w ij * s the total strengths of the 
network and 5 CiiCj is 1 if nodes % and j are in the same community and otherwise. 

In this paper communities are uncovered by optimising the modularity function in equation (2). 
The optimisation of Q is performed by using a tabu search algorithm [20] . We shall go back to some 
critical remarks on the use of modularity-based community-detection algorithms in the concluding 
Section. 

As discussed in Section I, one of the contribution of this paper is to compare commodity-specific 
community structures with a proper number of community benchmarks (as detailed in the next 
Section). To compare community partitions we use the normalised mutual information (NMI) 
measure, as introduced in [10]. To define the NMI index, we need to introduce the confusion 



matrix. Given two community partitions Va and Vb, the confusion matrix Af is defined as a 
matrix whose TVjj-th element is the number of nodes in the community % of the partition Va that 
appear in the community j of the partition Vb- The NMI is defined as: 



where Ca and Cb are respectively the number of communities in Va and Vb, N l = J2jNij, 
N.j = Nij and iV = £V ^ij- The NMI index is equal to 1 if Va and Vb are identical and 
assumes a value of if the two partitions are independent. 

IV. RESULTS 

A. Detecting the Community Structure of the Multi ITN 

We begin by studying the connectivity of the multi ITN and the size and concentration of its 
community structures. All results refer to the aggregate ITN and to the 14 commodity-specific 
layers as defined in Section II. In Table II we show the evolution of the density of the aggregate ITN 
(computed as the ratio between the number of existing links to the number of all possible links, 
i.e. N(N — 1)) and the relative density of commodity-specific networks (relative to the density 
of the aggregate ITN). We observe a monotonic increase in time of the aggregate ITN density, 
whereas the relative densities remained almost constant over time in each commodity-specific ITN. 
This implies an increase in the absolute value of commodity-specific densities. We also observe 
a relatively high heterogeneity of relative densities across commodity networks, which are always 
and significantly smaller than the aggregate one. This signals that results obtained using the 
aggregate ITN may be very different from those obtained looking at single commodity-specific 
networks (see also below). 

While the density measures the concentration of trade links in a network, the size of the largest 
connected component (LCC) measures its overall level of connectivity. In Table III, we report the 
size of the LCCs of the aggregate and commodity specific ITNs. While the former is always a 
completely connected network, this is not always the case for disaggregated cases, see for example 
arms and cereals. Other commodities, including electronics, optics, plastics and coffee, show 



NMI(Va,V b ) 




(3) 



7 



instead a large connectivity close to that of the aggregate ITN. This means that their contribution 
to overall connetivity is very strong. Notice also that for all the commodities we observe an 
increase in time in the size of the LCC, which is clearly a sign of the increase in the degree of 
integration of world trade. The largest changes in the size of LCC are observed for arms (c = 92), 
cereals (c = 10) and pharmaceutical products (c = 30). 

We now detect the community structure of both aggregate and commodity-specific ITNs by 
maximizing a weighted-directed version of the modularity function (see Eq. 2). We employ the 
community structure of the 12 yearly aggregate trade networks with weight matrix as in equation 
(1) as a first benchmark, in order to compare commodity-specific clusters with that obtained from 
the aggregate trade flows. 

The number of communities that we identify in each year and network is shown in Table 
IV. To begin, notice how the aggregate ITN typically displays a smaller number of communities 
than most of commodity-specific networks, meaning that the latter are more fragmented as far 
as trade clusters are concerned. In addition, we also observe that the number of communities in 
the aggregate ITN steadily increases over time, whereas this is not the case for most commodity- 
specific trade networks. In general, it appears that the smaller the size of the LCC, the higher 
the number of communities one finds. However, if for every ITN we look at the correlation across 
time between size of LCC and number of communities, results are different. While for some 
commodities a larger LCC size implies less communities, for others the opposite holds. In the 
first group we find coffee and tea (c = 9), pharmaceutical products (c = 30), precious stones 
(c = 71), and electric machinery (c = 84). In the second group we find all other commodities and 
the aggregate ITN. This evidence points to the existence of a large degree of heterogeneity in the 
number of community structures across commodity-specific networks, suggesting that the results 
obtained in the case of the aggregate ITN hide a lot of variability in the community structure of 
commodity-specific networks. This result is in line with similar one obtained in Ref. [7], where it 
is shown that many properties of the aggregate ITN (e.g., log- normal distributions of link weights 
and node-specific characteristics like strength and clustering) are the sheer result of aggregating 
their counterparts across heterogeneous commodity-specific networks. 

We now turn to a more detailed analysis of the community structure at a commodity-specific 
level. Figure 1 shows the distributions of the cluster size in year 2003. Again, the shape of the 
distributions and their ranges vary a lot across commodities. The commodities that generate 
the most concentrated community structures are electric machinery (c = 84), optical instruments 
(c = 89), and vehicles (c = 86), i.e. products that require more scientific knowledge. To draw a 
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more quantitative implication linking products and concentration of cluster-size distributions, we 
compute the normalized Herfindahl index (H), a synthetic measure of concentration of cluster size 
distributions. The index H, for a given commodity c and in a given year t, is defined as: 



H t C = T 



N 



N 



1 

N 



t = 1992,..., 2003, (4) 



where nx^,t is the number of communities identified in the network X t c and m c t {i) is the number 
of countries in the i-th. community in year t for commodity c. The index ranges between (no 
concentration at all) and 1 (maximum concentration). Table V reports the values of H£ for all 
networks and time periods. It is easy to notice that for the aggregate ITN there has been a decrease 
in concentration over time. This may be interpreted as a sign of the globalization process, as this 
pattern suggests that an increasing number of countries are participating to world trade over time. 
Indeed, while in 1992 we observe only 2 communities of about 80 countries each (one with Europe, 
Russia and Africa, the other with America and Asia), in 2003 a new community emerges, driven by 
China and India. At the commodity-specific level, an increase in H is observed also for coffee and 
tea (c = 9), mineral fuels (c = 27), pharmaceutical products (c = 30), arms (c = 92). However, for 
some other commodities we observe a decrease in H over time, see e.g. organic chemicals (c = 29), 
plastics (c = 39), and cotton (c = 52). This means that trade for those commodities has become 
less and less centralized and increasingly occurred among smaller and more dispersed groups of 
countries. 



B. Describing Trade Communities 

A useful way to visually describe community structure in the ITN is to employ colored world 
maps, where countries belonging to the same communities are associated to the same color. Figures 
2-4 report world maps depicting the community structure detected in 2003, for both the aggregate 
ITN and for the 14 commodity-specific networks. 

Notice that this visual device also allows us to informally correlate community structures with 
geographical considerations (we shall go back to a more formal analysis of this issue below). For 
example, most of the networks studied exhibit the presence of an American cluster composed of 
US and Canada (and often linked to Latin America), a European cluster (sometimes connected 
to North Africa), an Asian cluster consisting of China (and in many cases of India, Indochina and 
Australia) and finally a Russian community (sometimes linked to the European cluster). Africa 
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and Middle East are often split, independently of the commodity examined, among the other 
groups. This already suggests that geographical (and socio-political) factors are very important 
to explain the formation of community structure in the ITN. 

Apart from the regularities above, commodity-specific community structures often differ in a 
relevant way among each other. In what follows, we highlight some economically-relevant features 
of aggregate and commodity-specific community structures in 2003. We focus on 7 commodity 
classes, those exhbiting the most economically relevant patterns (the remaining 7 classes did 
not show such explicit regularities). Due to the relatively strong persistence over time of ITNs 
topological architecture (see Refs. [13] and [7] for a discussion), similar considerations also hold 
for other years. 

1. Aggregate ITNs: The world is divided in three major communities which follow a geo- 
graphical pattern: i) North and Latin America, ii) Europe, Russia, and North Africa, iii) 
China, India, Japan, Middle East, Australia and Sub-Saharan Africa. Two exceptions con- 
cern Africa: Nigeria and Ghana belong to the American community and we observe a minor 
separated community containing Belgium and the Democratic Republic of Congo, a former 
Belgian colony. 

2. Coffee and tea: We identify two communities containing coffee drinking and producing 
countries: i) Europe, Brazil, Peru, and Central African countries, ii) North America, Central 
America, Colombia and Venezuela. We also identify two communities of mainly tea drinking 
and producing countries: i) United Kingdom, South African and North-East African coun- 
tries, Pakistan and Bangladesh. Finally, there exist two mixed communities, but probably 
more connected with tea trade: i) India, Middle East, Russia, Australia, Argentina, Chile, 
ii) China, Japan, and Indochina. 

3. Cereals: The big producers of cereals belong each to a separate community: i) North 
America, ii) South America, iii) Russia. China and India are in separate communities, 
too. Finally, it is interesting to notice that Europe belongs to yet another separate cluster. 
Despite being not a big producer, but a big consumer, Europe is not an open market for 
agricultural products. This finding may be linked to the protectionist agricultural policies 
of the European Community. 

4. Mineral fuels: China and India have tight links with Middle East, Europe has links with 
Russia and North Africa, Brazil with Nigeria, and North America with Norway, which is 
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one of the largest oil producers in the world. 

5. Precious stones: In this case America and China belong to the same community. Europe, 
Russia, North and South Africa are the members of the largest community. Interestingly, 
countries rich of diamonds as Democratic Republic of Congo, Angola, and Sierra Leone, 
belong to a unique community containing also Israel. Finally, Australia, Indonesia, and 
India belong to another cluster. 

6. Electric machinery: There are only two communities strictly related to geographical dis- 
tance. Indeed, countries within a community share common borders. One contains America, 
China, Japan, India and Australia. The other one contains Europe, Russia, Africa and the 
Middle East. 

7. Vehicles: The world market for vehicles has a huge community containing America, China, 
Japan, India, Australia and almost all African countries. This may reflect the high diffusion 
of Japanese cars in Africa. Russia is still a closed market containing all former Soviet 
republics. Finally, Europe is divided in two communities, which have almost no members 
outside the continent: a finding that seems to reveal a protectionist market for vehicles in 
Europe. 

8. Arms: The community structure for arms is highly fragmented and therefore difficult to 
interpret. Moreover, many countries seem not to belong to any community. Interestingly, 
these countries are often those were civil wars or in general social instability are most likely 
to be (or to have been) present. This is the case of Mozambique, Zambia, Angola, Guinea, 
Myanmar and Central Asian countries. It is unlikely that these countries do not participate 
in arms trade, but it is not surprising that our data do not reveal this, as probably that kind 
of trade relationships are not official. Finally, Africa is the most fragmented continent and 
almost all communities found contain some African countries. 

C. Comparing Community Structures 

In this Section, we explore more quantitatively commodity-specific community structures by 
using the NMI index introduced in Section III. 

To begin with, we ask to what extent community structures are stable over the time interval con- 
sidered. To do that, we compare the partitions obtained at time t and t + 1 for t — 1992, . . . , 2002. 
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More precisely, for the aggregate ITN and for any c, we compute the quantity NM1(V£,V£ +1 ), 
where is the partition of our iV countries in year t for commodity c. This gives a measure of 
stability over time of community structures (see Table VI). Notice how the smallest values of NMI 
(large community structure changes) are observed in the early 1990s. In more recent years, on the 
contrary, NMIs have been larger, meaning weaker changes in the composition of communities from 
year to year. If one instead compares partitions in 1992 with those in 2003 (i.e., one computes the 
quantity NMI(7 :, f 992 , P| 003 )), it turns out that the stronger changes are associated to coffee and 
tea (c = 9), pharmaceutical products (c = 30), and arms (c = 92). The most stable community 
structures are instead those of aggregate trade, plastics (c = 39), optical instruments (c = 89), 
mineral fuels (c = 27), iron and steel (c = 72), and cotton (c = 52). Notice also that, on aver- 
age, the majority of commodity-specific community structures were less stable than that of the 
aggregate network. Again, this suggests a strong mismatch between aggregate and disaggregated 
properties. 

We then compare, in each given year, partitions associated to the aggregate ITN with those 
associated to commodity-specific networks by computing the quantities NML(Vf l ,V£), where Vf l 
is the partition obtained from the aggregate ITN. This exercise is meant to ask more quantitatively 
the question whether the aggregate community structure can well predict those obtained at the 
commodity-specific level, or, put it differently, the extent to which the community structure of 
any commodity-specific network contributes to shape (or is able to predict) that observed at the 
aggregate-trade level. Inspection of Table VII shows that NMI values are increasing in time for 
almost all commodities. This means that commodity-specific community structures are becoming 
more and more similar to the aggregate one, i.e. that the role of all commodities in shaping 
the aggregated community structure has increased in time. In particular, we observe the largest 
increase in NMI from 1992 to 2003 for mineral fuels (c = 27), plastics (c = 39), iron and steel 
(c = 72), pharmaceutical products (c = 30), and organic chemicals (c = 29). Moreover, for all 
the years considered, mineral fuels and plastics appear to be the commodities whose community 
structure is the most similar to the aggregated one. Overall, a major role for the chemical sector 
emerges from these results, in that they are the ones whose country partition better mimics the 
aggregate one. 

The overall similarity pattern between commodity-specific partitions can also be studied using 
minimum spanning tree (MST) techniques. Indeed, starting from the similarity between commod- 
ity i and commodity j expressed by the index NMI(i, j), we can define a distance between two 
community structures as 1 — NMI(i, j). Using this metrics, we can build a MST as in Ref. [24] 
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in order to classify the commodities into groups displaying the highest similarity in terms of their 
community structure. The results for year 2003 is in Figure 5 (similar results hold for all the 
remaining years). The figure shows that commodities related to science- or technology-based in- 
dustries (nuclear reactors, optical instruments, electric machinery) are the most similar in terms of 
their community structures, whereas arms appears to be the most dissimilar one. More generally, 
this exercise shows that it is possible to find meaningful classifications of commodities using in- 
dicators assessing the similarity between community structures characterizing commodity-specific 
trade networks. 

D. Community Structure, Geography, and Trade Agreements 

We finally turn to study the extent to which community structures identified using trade data 
correlate with other economics-relevant data. To do that, we employ the geographical closeness 
matrix S, as well the time-dependent RTA networks M t . The entries of the symmetric and time- 
independent matrix S, to repeat, express a measure of geographical closeness between pair of 
countries, computed as the inverse of geographical distance between their most populated cities. 
The entries of the symmetric but time-dependent matrices M t contain, in a give year, the number 
of trade agreements currently in place between any two countries, irrespective of the type of RTA 
signed (bilateral, multilateral, commodity- specific, etc.). The underlying assumption is that the 
higher this number, the closer the two countries are in the RTA space (and thus, according to 
empirical findings, the larger their expected trade flows). 

We apply to both S and M t the community-detection algorithms explained in Section III 
and previously applied to trade matrices. Therefore, we end up with a geographically-induced 
community partition J> GEO and 12 time-dependent RTA-induced partitions Vf' . The resulting 
partitions are visualized in the maps of Figure 6. Notice that clusters in p GEO represent groups 
of countries that are geographically close, without using exogenously-determined partitions of 
countries (e.g., based on continents or sub-continental breakdowns). The community structures 
in J>^ TA do instead pick up clusters of countries that not only belong to free-trade or multilateral 
agreements (e.g. NAFTA, Mercosur, EU, etc.), but also signed additional bilateral agreements. 

We first compare, using the NMI, the aggregate ITN community structure with those detected 
using geographical distances or RTAs (see Figure 7). We observe increasing NMIs across time 
until 2001 and a slight decrease afterwards. We also find more similarity between aggregate 
trade and geography based communities with respect to communities determined by RTAs. Thus, 
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geographically-related factors seem to explain the pattern of global trade more than political de- 
terminants. Also, this result is more evident in the recent years after 2001. A possible explanation 
might be the global political crisis after 11th September 2001 that implied a slight decrease in 
global trade as a consequence of the wars in Iraq and Afghanistan. 

When comparing community structures of commodity-specific ITNs with the partitions ob- 
tained from geographical and RTA data (see Tables VIII and IX), we find results similar to the 
aggregate case. In general, it is geography and not trade agreements that seems to correlate more 
with the observed patterns. Plastics (c = 39) and mineral fuels (c = 27) display the highest 
similarity with RTA communities. The same result holds when confronting trade communities 
with geographical data, but in addition we notice high NMIs also for iron and steel (c = 72) and 
cotton (c = 52). 

These results reinforce the traditional view put forth by standard gravity-equation trade em- 
pirics [11], which stresses the importance of geographical distance (as a proxy for trade resistance 
factors) in determining bilateral trade flows. Here, we show that geographical distance is impor- 
tant to predict not only the expected flow of a bilateral trade relationship (e.g., exports from 
country A to country B), but also the formation of trade communities, that is complicated trade 
structures multilaterally involving groups of countries. On the other hand, our findings contribute 
to the discussion related to the impact of international agreements on world trade and seem to go 
in the direction of Ref. [32], which shows that there is no evidence that the WTO has increased 
international trade. 

V. CONCLUDING REMARKS 

In this paper, we provide a first exploratory study of the community structure of commodity- 
specific trade networks from 1992-2003. After recovering the optimal partition of countries, we 
compare commodity-specific communities with the aggregate trade community. 

Our results show that commodity-specific community structures are very heterogeneous and 
in general their statistical properties are quite different from those of the community structure 
of the aggregate ITN. For example, whereas the number of communities of the aggregate ITN 
increases in time, this is not the case for most commodity-specific trade networks. Moreover, the 
shapes and ranges of cluster-size distributions vary a lot across commodities. As far as community 
structure evolution is concerned, one observes a decrease in concentration over time of cluster-size 
distributions (a sign of the globalization process) for the aggregate ITN, a pattern that is not always 
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matched at the commodity-specific level, where trade associated to some products have become 
less and less centralized and increasingly occurred among smaller and more dispersed groups of 
countries. Furthermore, the community structure of the aggregate ITN has been changing more 
slowly over time than their commodity-specific counterparts. 

We have also explored to what extent the community structure of any commodity-specific 
network may contribute to shape (or is able to predict) that observed at the aggregate-trade 
level. We have shown that commodity-specific community structures are becoming more and 
more similar to the aggregate one, i.e. that the role of all commodities in shaping the aggregated 
community structure has increased in time. However, a major role for the chemical sector appears 
from these results, in that they are the ones whose country partition better mimics the aggregate 
one. More generally, by comparing the commodity-specific community structures using the NMI 
indicator, one can find meaningful classifications of commodities. 

Finally, we have explored two possible factors that correlate with community structure, namely 
geographical distance and the existence of regional trade agreements between countries. Our 
findings suggest that geographical distance correlates much more with the observed community 
structure than RTAs. This result confirms previous findings from the empirical literature on trade. 

The paper can be extended and refined in at least three directions. First, our findings related 
to the impact of geography and RTAs are only partial, as they only check for unconditional effects 
(i.e. they do not address the residual effects of trade agreements once geography is controlled for). 
In order to make our statements more robust, one may follow Ref. [31] and compare communities 
observed in trade data with those detected in the network built with the predictions of a standard 
gravity model [see for example Ref. 11] . 

Second, the robustness of our results should be checked against a number of possible problems. 
For instance, it is well-known that modularity-based community detection suffers from a resolution 
limit bias [16]. Therefore, community detection algorithms based on alternative criteria may 
be employed (e.g., community detection methods based on information theory, see for example 
Ref. [33]). Similarly, one may consider to apply algorithms allowing for overlapping communities 
[27, 28]. 

Third, another point that deserves further analysis is the detection of community structures 
across commodity-specific layers. In the paper, we have analyzed independently the most impor- 
tant 14 layers. This allows one to identify groups of countries that trade the same commodity 
among them. From an economic point of view this signals strong interdependencies but does not 
convey any insights on the input-output structure of the cluster. For example, there might be 
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groups of countries that are linked in tightly connected chains or cycles, where a country imports 
from another one a particular type of commodity needed as input for its peculiar industrial struc- 
ture, and at the same time exports to other countries in the group another commodity that is 
fed into their production processes (or consumed as final good). In order to address these issues, 
one would like to either synthesize into a meaningful statistic all commodity-specific relationships 
between any two countries or apply new techniques able to detect community structures in multi 
graphs [25]. 



[1] The choice of a 2-digit breakdown of data may be considered insufficient to clearly identify homo- 
geneous product lines, but it has been made because in the HS classification system there is not 
a unique way to further disaggregate flows by commodities at a higher number of digits. Notice, 
however, that network analyses often face a trade off between the need for a finer disaggregation and 
the very possibility to obtain connected graphs: typically, as soon as 3 or 4 digit data are considered, 
the resulting graphs easily become not connected, with the size of the largest connected component 
quickly decreasing. 

[2] Since, as always happens in trade data, exports from country i to country j are reported twice 
(according to the reporting country — importer or exporter) and sometimes the two figures do not 
match, we follow Ref. [14] and only employ import flows. For the sake of exposition, however, we 
follow the flow of goods and we treat imports from j to i as exports from i to j. 

[3] We refer the reader to Ref. [7] for a thorough analysis of the topological properties of this database 
from a multi-network perspective. 

[4] See http://www.wto.org/. 

[5] Results are robust to alternative distance measures. Data are available at the URL: 
http://www.cepii.fr/. 

[6] A. Arenas, J. Duch, A. Fernandez, and S. Gomez. Size reduction of complex networks preserving 
modularity. New Journal of Physics, 9(6):176, 2007. 

[7] M. Barigozzi, G. Fagiolo, and D. Garlaschelli. Multi-network of international trade: A commodity- 
specific analysis. Physical Review E, 81:046104, 2010. 
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Code 


Commodity 


Value 
(USD) 


Value 
(USD) 


per Link 


% of Aggregate 
Trade 


83 


Nuclear reactors, boilers, machinery and 
mechanical appliances; parts thereof 


5.67 x 10 11 


6.17 x 


10 7 


11.37% 


84 


Electric machinery, equipment and parts; 
sound equipment; television equipment 


5.58 x 10 11 


6.37 x 


10 7 


11.18% 


27 


Mineral fuels, mineral oils & products of their 
distillation; bitumin substances; mineral wax 


4.45 x 10 11 


9.91 x 


10 7 


8.92% 


86 


Vehicles (not railway, tramway, rolling stock); 
parts and accessories 


3.09 x 10 11 


4.76 x 


10 7 


6.19% 


89 


Optical, photographic, cinematographic, 
measuring, checking, precision, medical or 
surgical instruments/apparatus; 
parts & accessories 


1.78 x 10 11 


2.48 x 


10 7 


3.58% 


39 


Plastics and articles thereof 


1.71 x 10 11 


2.33 x 


10 7 


3.44% 


29 


Organic chemicals 


1.67 x 10 11 


3.29 x 


10 7 


3.35% 


30 


Pharmaceutical products 


1.4 x 10 11 


2.59 x 


10 7 


2.81% 


72 


Iron and steel 


1.35 x 10 11 


2.77 x 


10 7 


2.70% 


71 


Pearls, precious stones, metals, coins, etc 


1.01 x 10 11 


2.41 x 


10 7 


2.02% 


10 


Cereals 


3.63 x 10 10 


1.28 x 


10 7 


0.73% 


52 


Cotton, including yarn and woven fabric thereof 


3.29 x 10 10 


6.96 x 


10 6 


0.66% 


09 


Coffee, tea, mate & spices 


1.28 x 10 10 


2.56 x 


10 6 


0.26% 


92 


Arms and ammunition, parts and 
accessories thereof 


4.31 x 10 9 


2.46 x 


10 6 


0.09% 


ALL 


Aggregate 


4.99 x 10 12 


3.54 x 


10 8 


100.00% 



TABLE I. The 14 most relevant commodity classes in year 2003 in terms of total-trade value (USD), 
trade value per link (USD), and share of world aggregate trade. 
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FIG. 1. Cluster-size distributions in 2003. 
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(a)Coffee and tea c = 9. (b)Cereals c = 10. 




(c)Mineral fuels c = 27. (d)Organic chemicals c = 29. 




(e)Pharmaceutical products c = 30. (f)Plastics c = 39. 




(g)Cotton c = 52. (li)Precious stones c = 71. 




(i)Iron and steel c = 72. (j)Nuclear reactors c = 83. 



FIG. 3. World maps showing trade communities of commodity specific ITNs in 2003. In gray countries 
not belonging to any community or for which no data are available. 



25 



(a)Electric machinery c = 84. 



(b)Vehicles c = 86. 




(c)Optical instruments c = 89. (d)Arms c = 92. 

FIG. 4. World maps showing trade communities of commodity specific ITNs in 2003. In gray countries 
not belonging to any community or for which no data are available. 
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FIG. 5. Minimum spanning tree for 2003. The minimum spanning tree (MST) is computed starting from 
the similarity between commodity % and commodity j expressed by the index NMI(i, j) and defining a 
distance between two community structures as 1 — NMI(z, j). Codes (from left to right): 83 = nuclear 
reactors; 89 = optical instruments; 84 = electric machinery; 39 = plastics; 52 = cotton; 72 = iron and 
steel; 27 = mineral fuels; 30 = pharmaceutical products; 10 = cereals; 29 = organic chemicals; 86 = 
vehicles; 71 = precious stones; 09 = coffee and tea; 92 = arms. 
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(a)RTAs. (b)Distances. 
FIG. 6. World maps showing RTAs in 2003 and geographic communities. 
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FIG. 7. NMI when comparing the community structures induced by the exogenous networks build using 
geographical distances (GEO) or regional trade agreements data (RTA) with the the community structures 
of aggregate trade. 
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